A block Chebyshev-Davidson method with inner-outer restart for large eigenvalue problems
نویسنده
چکیده
We propose a block Davidson-type subspace iteration using Chebyshev polynomial filters for large symmetric/hermitian eigenvalue problem. The method consists of three essential components. The first is an adaptive procedure for constructing efficient block Chebyshev polynomial filters; the second is an inner–outer restart technique inside a Chebyshev–Davidson iteration that reduces the computational costs related to using a large dimension subspace; and the third is a progressive filtering technique, which can fully employ a large number of good initial vectors if they are available, without using a large block size. Numerical experiments on several Hamiltonian matrices from density functional theory calculations show the efficiency and robustness of the proposed method. 2010 Elsevier Inc. All rights reserved.
منابع مشابه
A block Chebyshev-Davidson method for linear response eigenvalue problems
We present a Chebyshev-Davidson method to compute a few smallest positive eigenvalues and corresponding eigenvectors of the linear response eigenvalue problem. The method is actually applicable to the slightly more general linear response eigenvalue problem where purely imaginary eigenvalues may occur. For the Chebyshev filter, a tight upper bound is obtained by a computable bound estimator con...
متن کاملRestarting iterative projection methods for Hermitian nonlinear eigenvalue problems with minmax property
In this work we present a new restart technique for iterative projection methods for nonlinear eigenvalue problems admitting minmax characterization of their eigenvalues. Our technique makes use of the minmax induced local enumeration of the eigenvalues in the inner iteration. In contrast to global numbering which requires including all the previously computed eigenvectors in the search subspac...
متن کاملRayleigh Quotient Iteration and Simplified Jacobi-davidson with Preconditioned Iterative Solves for Generalised Eigenvalue Problems
The computation of a right eigenvector and corresponding finite eigenvalue of a large sparse generalised eigenproblem Ax = λMx using preconditioned Rayleigh quotient iteration and the simplified JacobiDavidson method is considered. Both methods are inner-outer iterative methods and we consider GMRES and FOM as iterative algorithms for the (inexact) solution of the inner systems that arise. The ...
متن کاملSome new restart vectors for explicitly restarted Arnoldi method
The explicitly restarted Arnoldi method (ERAM) can be used to find some eigenvalues of large and sparse matrices. However, it has been shown that even this method may fail to converge. In this paper, we present two new methods to accelerate the convergence of ERAM algorithm. In these methods, we apply two strategies for the updated initial vector in each restart cycles. The implementation of th...
متن کاملA Parallel Scalable PETSc-Based Jacobi-Davidson Polynomial Eigensolver with Application in Quantum Dot Simulation
The Jacobi-Davidson (JD) algorithm recently has gained popularity for finding a few selected interior eigenvalues of large sparse polynomial eigenvalue problems, which commonly appear in many computational science and engineering PDE based applications. As other inner–outer algorithms like Newton type method, the bottleneck of the JD algorithm is to solve approximately the inner correction equa...
متن کاملذخیره در منابع من
با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید
برای دانلود متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید
ثبت ناماگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید
ورودعنوان ژورنال:
- J. Comput. Physics
دوره 229 شماره
صفحات -
تاریخ انتشار 2010